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Abstract. Arguments of astronomers against equation of motion for arbitrarily shaped 
particle under action of electromagnetic radiation are discussed. Each of the arguments is 
commented in detail from the point of view of the required physics. It is shown that the 
arguments of astronomers, including referees in several astronomical and astrophysical 
journals, are unacceptable from the physical point of view. Detail explanations should help 
astronomers in better physical understanding of the equation of motion. Relativistically 
covariant equation of motion for real dust particle under the action of electromagnetic 
radiation is derived. The particle is neutral in charge. Equation of motion is expressed in 



X 

terms of particle's optical properties, standardly used in optics for stationary particles 
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1. Introduction 

Astronomers each year publish several papers on orbital evolution of cosmic (dust) par- 
ticles which consider interaction of the particles with electromagnetic radiation. This is 
taken in the form of the Poynting- Robertson effect (P-R effect; Robertson 1937). How- 
ever, P-R effect corresponds to a very special form of interaction between the particle 
electromagnetic radiation (see Eqs. (120), (122) in Klacka 1992a). Since real particles do 
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not fulfill this special form of interaction, one needs to have more general equation of 
motion in disposal. This is presented in Klacka (2000a). Later on, more simple deriva- 
tions were presented: Klacka (2000b), Klacka and Kocifaj (2001a), sec also review paper 
Klacka (2001b). Applications of the equation of motion can be foimd in Klacka (2000c), 
Klacka and Kocifaj (2001a), Klacka and Kocifaj (2001b), Klacka and Kocifaj (2002), 
Kocifaj and Klacka (2002a, 2002b). 

Our experience shows that astronomers do not understand the physics of the equation 
of motion. This is the reason why we have decided to present arguments against the 
equation of motion and, of course, to give detail explanations why the arguments are 
physically incorrect. We can say that this paper is a continuation of the paper Klacka 
(1993). 

2. Arguments and answers 

This section presents arguments/statements of astronomers against the equation of mo- 
tion (since the year 1998) and our point of view. 

2.1. Argument/statement 1 

Astronomers do not need any more general equation of motion - they are satisfied with 
the Poynting-Robertson effect. 
Answer: 

Scattering of light, electromagnetic radiation, on arbitrarily shaped (dust) particles 
with various optical properties may significantly differ from that required by the P-R 
effect: compare Klacka (1992a) and Kocifaj and Klacka (1999). As a consequence, orbital 
evolution of real particle practically always differs from that corresponding to the P-R 
effect: Kocifaj et al. (2000), Klacka and Kocifaj (2001a), Klacka and Kocifaj (2001b), 
Klacka and Kocifaj (2002), Kocifaj and Klacka (2002a, 2002b). 

2.2. Argument/statement 2 

Orthonormality of unit vectors used in the scattering theory holds also in the frame of 
reference in which particle moves with velocity v ( also for v 
Answer: 

Although this access was used in Klacka (1994), Klacka and Kocifaj (1994) and Kocifaj 
et al. (2000), it is not physically correct. Unit vectors must be orthonormal in the proper 
frame of reference (the rest frame of the particle, i. e. a moving frame relative to the 
source of electromagnetic radiation). The requirement that "final equation of motion has 
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to be written in a covariant form" offers only one possibility: vectors are not orthonormal 
in any other frame of reference. And physical interpretation is also evident - aberration 

of light (Klacka 2000a, 2000b). Thus, the statement about the conservation of normality 
of vectors under Lorentz transformation is nonphysical. 

However, experience shows that astronomers do not accept this general argument. 
Thus, we will present a simple example which shows that access of Klacka (2000a) is 
correct (at least for the simple example). 



EXAMPLE: 



Let us consider a plane mirror moving (at a given moment) along x-axis (system S) 
with velocity v = {v,0,0), v > 0; the mirror is perpendicular to the x-axis (the plane of 
the mirror is parallel to the yz-plane). A beam of incident (hitting) photons is character- 
ized by unit vector S' = {cos6' ,sm9' ,0) in the proper frame (primed quantities) of the 
mirror. Reflected beam is described by the unit vector e' = (— cos sin ^',0) (in the 
proper frame S'). 



The problem is: Find equation of motion of the mirror in the frame of reference S. 



SOLUTION 1: trivial manner 



Consider one photon (frequency /') in the proper frame of the mirror. Since the 
directions (and orientations) of the incident and outgoing photons are characterized by 

S' = {+ cos 6*', sin 6*', 0) , 

e' = (- cos 6'', sin 61', 0) , (1) 

we can immediately write 

p'f = (1,+ cos^',sin0',O) , 

h f 

= -i- (1,- cos 6*', sin 6*', 0) , (2) 

for the four-momentum of the photon before interaction with the mirror and after the 
interaction. 

As a consequence, the mirror obtains four-momentum 

P^=pt - p'o^ = — (0,2 cos^',0,0) . (3) 



Application of the special Lorentz transformation to Eq. (3) yields 
p'* = ^27 (cos^O (A1,0,0) , (4) 

where, as standardly abbreviated, 
7 = 1/ sjl - /32 ■ (i = v I c. (5) 

On the basis of Eq. (4), we can immediately write equation of motion of the mirror 

^ = ^27 (cos e') (A 1,0,0) , (6) 
At c 

where E[ is the total energy (per unit time) of the incident radiation measured in the 
proper frame of reference. 

SOLUTION 2: application of general theory of Klacka (2000a) 

We have to choose orthonormal vectors in the systems S': we will use S' and e'^ and 
one can easily find 

S' = (+ cos 6*', sin 6*', 0) , 

e'l = (- sin 61', cos 6*', 0) . (7) 
We have to write {Q'2 = 0) 

P'=^- (Q'r + Qi e'l) • (8) 
On the basis of Eqs. (3), (7) and (8) we have 

= 2 (cos6'')^ , g'l = - 2 (sin6»') (008 6*') ; (g^ = 0) . (9) 
Other prescription yields 
6° = 7 (1 + . S'/c) = 7(1 + /? cos^O , 

6- = S' + [(7 - l)v-S'/v'^ + 7/c] t; = (7 cos 6'' + 7 /3,sin6»',0) , (10) 
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= j{l + V e[/c) =7(1-/3 sin 61') , 

= e'l + [(7 - 1) V -e'^/v^ + 7/c] v = {- J sine' + 7 /3,cos6'',0) 



(11) 
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Inserting Eqs. (9) - (11) into Eq. (28) of Klacka (2000a), one obtains 
Hn^^ E'. 

^ = {[2(cos^') ] {b'^ - + [-2(sin^') (cos^')] (&i - P^} 
E' 

= -^27 (cos6'') (/3, 1,0,0) . (12) 



COMPARISON: 

Unit vectors <S" and e'l are used according to Klacka (2000a). They lead to correct 
results. The unit vectors are orthonormal in the system S'. How does the situation look 
in the system S? 

S = ^ {S' + [{j ~ l)v S'/v' + 7/c] v} , 
w 

w' = ^ {1 + vS' / c) , (13) 

and analogous equation holds for vector ei. Inserting Eqs. (7), one obtains 

_ cosg' + /? sin6>^ 
^ ^ I + P cose*'' 7(1 + /3 COS0')' 



r - sin6l' + j3 ^ , 

"^ = | l-/3sin^^ ' 7(l — ^''^ 



cose' 



It can be easily verified that scalar product of these two vectors is nonzero, in general, 
even to the first order in /?: 

S-ei «/? (cos6'' - sin 6*') . (15) 



SOLUTION 3: application of general theory of " bi ^ (Klacka 2001a). 
We will use equation of motion in the form: 

f = ^Q'S'. (16) 
dr c 

Calculations yield 

Q' = Qi - Q'r= diag{l, 1, 1) - diag{- 1, 1, 1) = diag{2, 0, 0) . (17) 
Moreover, expression for S presented in Eq. (14) yields 

\1 + p cosO' ) 

{Q'&f = 2 (/?,0,0) , (18) 
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0^{Q'S) = 2/3 (cos 6*' + /3) / (1 + /? cos6i') 
/3^(g'/3) = 2 0" . 

Inserting into equation 

, — ^ V 1 

(XT C 

one obtains 
^ = ^27 (cos^O (/?, 1,0,0) . 



(19) 
(20) 
(21) 



or c 

REMARK: Inserting i;,' = S A' cos^' into Eq. (6) (or Eqs. (12), (21)) and using Eq. 
(14) for the purpose of obtaining cos^' = (cos^ — /3)/(l — /? cos^), one easily obtains: 

i) dE/dT = 2^^ S A' {cos 9 — (3)^ /3; using definition of radiation pressure dE/dt = 
P V A', we have P = 2 {S / c) (cos 61 - f3)^ / (1 - /J^), or, 

ii) dp/dr = 2 {S A' / c) (cos 6 — f3}^; using definition of radiation pressure P = 
{dp/dt) I A', wc have P = 2 [S / c) {cos6 - (if / (1 - 13^). 

Result for P is consistent with result presented in Einstein (1905). 

2.3. Argument/ statement 3 

The assumption is that energy E' of the particle is unchanged: the energy of the incoming 
radiation equals the energy of the outgoing radiation, per unit time. What this means 
is that the particle does not heat up. Fine, as such, but wait. A physical explanation 
here involving the Poynting-Robertson effect has the particle losing orbital energy due to 
re-radiation effects, and so the size of its orbit reduces. So how does it work, physically? 
Answer: 

Yes, dE' / dr = corresponds to conservation of mass of the particle. However, the 
relation for energy holds in the proper reference frame of the particle (rest frame of 

the particle), only. Lorcntz transformation yields that energy changes in other reference 
frames, e. g. in the rest frame of the source of radiation. The special case Q[ = Q'2 = 
in Klacka (2000a - c. g., Eqs. (7), (23)) - corresponding to the P-R effect yields 
that dE/dt 7^ and the corresponding change of semi-major axis is given by Eq. (14) in 
Klacka (1992b). 
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2.4- Argument/ statement 4 

Where does the "really [realistically?] shaped particle come into it?" 
Answer: 

Effective factors Q'^, Q[, QJ, (see Eq. (7) in Klacka 2000a) take into account optical 
characteristics of the particle. They can be calculated using, e. g., discrete dipole approx- 
imation - see, e. g., Draine and Weingartner (1996). 



2.5. Argument/statement 5 

dE' / dr = 0: If a momentum is transferred to the particle there should be a corresponding 
gain in kinetic energy. In reality this term will depend upon the particle's albedo etc. 

Answer: 

The statement dE' / dr = was discussed in section 2.3 - conservation of particle's 
mass. Even within a Newtonian physics we have T = m v'^/2, dT/dt = m v ■ v and this 
yields dT'/dt = in the proper frame of the particle, although dp' /dt ^ 0. In reality 
dE/dt ^ in the rest frame of the source and dE/dt depends on optical properties of the 
particle (see sections 2.3 and 2.4). (It has no sense to use the term "particle's albedo" if 
the size of the particle is comparable to the wavelength of the incident electromagnetic 
radiation, when measured in the proper frame of the particle.) 



2.6. Argument/ statement 6 

The problem is of great general interest and (at least at the Newtonian level) has been 
treated intensively in the literature. 

Answer: 

Since the time of Maxwell and Einstein (1905) we know that the discussed problem 
cannot be treated at the Newtonian level. In reality, motion of cosmic dust grains is stan- 
dardly described by the P-R effect (of course, it cannot be imderstood at the Newtonian 
level). P-R effect is a very special case of the general equation of motion (Eq. (24) in 
Klacka 2000a). 
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2.1. Argument/ statement 7 

A real astronomical problem is oversimplified and the resulting equations are not really 
useful. 

Answer: Papers Klacka and Kocifaj (2001a, 2001b) show important difference between 
orbital motion obtained on the basis of Eq. (24) in Klacka (2000a) and the P-R effect. 
P-R effect was used for six decades (up to now) and there was no astronomical protest 
that " it is not really useful" . 

3. Conclusion 

We have discussed some statements of astronomers as for more general equation of motion 
than that given by the P-R effect. We have tried to "stand up for" the general equation 
of motion, as we have presented it for the first time in Eq. (24) in Klacka (2000a) and 
in a more simple form in Klacka (2000a), in Klacka and Kocifaj (2001a). We have shown 
that our more general equation of motion reduces to the two cases discussed earlier in 
literature: Einstein (1905), Robertson (1937). Any correct equation of motion has to be 
consistent with the results of Einstein (1905) and Robertson (1937). 
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